Covariance for Conic and Wedge Complete Filling. 
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Interfacial phenomena associated with fluid adsorption in two dimensional systems has recently 
been shown to exhibit hidden symmetries, or covariances, which precisely relate local adsorption 
properties in different confining geometries. We show that covariance also occurs in three dimen- 
sional systems and is likely to be verifiable experimentally and in Ising model simulations studies. 
Specifically, we study complete wetting in wedge (W) and cone (C) geometries as bulk coexistence 
is approached and show that the equilibrium mid-point heights satisfy lc{h, a) — lw{^ , ct), where h 
measures the partial pressure and a is the tilt angle. This covariance is valid for both short-ranged 
and long-ranged intermolecular forces and identifies both leading and next-to-leading order critical 
exponents and amplitudes in the confining geometries. Connection with capillary condensation-like 
phenomena is also made. 
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The central result of density functional theory, which 
underpins the modern description of inhomogeneous flu- 
ids, is that the external potential is a unique functional 
of the equilibrium density profile p(r). That is, p{r) is 
sufficient to completely determine both the range of the 
surface intermolecular forces and the shape of the con- 
fining walls P, Q . 

In spite of this rigourous result, recent studies have 
highlighted precise connections, or covariances, between 
adsorption properties for fluids at differently shaped sub- 
strates. For instance, the adsorption occuring at wedges 
and apexes can be related to adsorption properties at a 
planar substrate in a non-trivial manner !3| . Covariance 
behaves like a hidden symmetry for fluid interfaces and 
has a number of profound implications for the theory of 
wetting, including the allowed values of critical exponents 
and the very structure of interfacial Hamiltonians 0]. 

In this paper, we describe a new example of geomet- 
rical covariance which relates the interfacial properties 
of fluids in three dimensional wedges and cones. Unlike 
previous cases, the covariance we describe here occurs 
off bulk coexistence, for arbitrary intermolecular forces 
and in the complete wetting regime, i.e. for zero con- 
tact angle. This is of particular relevance since complete 
wetting is readily accessible in the laboratory 5] , imply- 
ing that the existence of geometrical covariances for fluid 
interfaces could be verified experimentally for the first 
time. Indeed, complete wetting in wedges has already 
been studied in the laboratory with high accuracy 1^ . 

The main result of this paper is the covariance relation 



lc{h, a) ^ Iwi - , a), 



(1) 



where lw{h,a) and lc{h,a) denote the equilibrium inter- 
facial heights at the mid-point of a wedge and a cone, 
respectively. Here a represents the substrate tilt angle 
whilst h measures the partial pressure difference Psat — P- 



This relation is valid as /i and describes all diverg- 
ing contributions to the film thicknesses. We emphasise 
that similar to other cases of covariance, the above re- 
lation exists despite significant differences in the phase 
transitions occuring in the two geometries. 

To begin, we review briefly the basic theory of com- 
plete wetting 0. Consider a planar wall in contact with 
a bulk vapour at temperature T and chemical potential 
fi. We suppose that T is above the wetting temperature 
so that as /i is increased towards saturation fisat{T), 
the thickness of the adsorbed liquid layer l-^ diverges, 
corresponding to zero contact angle 9 — 0. The diver- 
gence is characterised by the power-law lTr{h) ~ h~^= 
and is accompanied by the divergence of the transverse 
correlation length ^|| ~ and interfacial roughness 

The standard way of modelling fluctuation effects at 
complete wetting is via the effective Hamiltonian Q 



Wil) 



(2) 



where E is the surface tension of the liquid-vapour inter- 
face. For systems with long-ranged forces, the binding 
potential W{1) is given by 



W{l) = hl 



A 



(3) 



where A is a Hamaker constant and p depends on the 
range of the forces. For the experimentally relevant case 
of three dimensions, the cases p = 2, 3 correspond to van 
der Waals (vdW) and retarded vdW forces respectively. 
In three dimensions, ffuctuation effects are not important 
and a mean-field analysis suffices to determine the critical 
properties yielding /?™ = l/{p+l), v\\ = {p -\- 2) / 2{p + 1) 
and ~ Y^ln^ii . In two dimensions, mean-field theory 
breaks down for p > 2 and the exponents are universal 
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/3™ = = 1/3 and = 2/3 corresponding to the so- 
called weak fluctuation regime Q- 

Now consider the analogous complete wetting in a 
three dimensional wedge. The wedge is characterised by 
a tilt angle a so that the height of the substrate above 
the horizontal (say) is ip{x,y) = tana|a;| (see Fig. 1). 
The wedge geometry enhances the adsorption so that 
lw{h, a) is far greater than the thickness of the layer ad- 
sorbed at a flat substrate l-yrih). The divergence of 1^ 
as /i ^ is accompanied by the divergence of corre- 
lation lengths £,y across and along the wedge, and 
also by the (mid-point) roughness ^j^. However, fluctu- 
ation effects are not particularly important at complete 
wedge wetting and all these length scales can be related 
to Iw In particular, ^ S.y ^ Iw and ^ y/lnlw The 
dominant divergence of Z^, follows from simple thermo- 
dynamic considerations since the macroscopic meniscus 
must be an inscribed cylinder with radius as de- 

termined by the Laplace pressure |^. Thus, at leading 
order, 1^ « H[seca~ l)/h. Similar critical properties oc- 
cur for the two dimensional wedge although for this case 

To go beyond the macroscopics and obtain a more ac- 
curate expression for 1^ , we resort to an effective interfa- 
cial Hamiltonian description This reveals the pres- 
ence of a significant next-to-leading-order contribution in 
Iw both in two and three dimensions which itself becomes 
macroscopic as /i — )■ 0. For convenience, we only present 
details based on the shallow wedge model, 

H[l,^]^ jd^i^^{Vlf+W{l-^)Y (4) 

appropriate when tana « a. Nevertheless, analysis of 
a more involved drum-head model results in identical fi- 
nal expressions which are valid for arbitrary tilt angles. 
Similar to complete wetting at a flat wall, a mean-field 
analysis suffices to determine the critical properties in 
three dimensions since both the wedge thickness 1^ and 
the wetting layer thickness Z^r are much greater than the 
roughness. The Euler-Lagrange equation for the relative 
height ri{x) = l{x) — a\x\ obtained from minimising the 
Hamiltonian is 

^V"^W'{^), (5) 

which is solved subject to the boundary conditions 
77'(0~'') = —a and rj'{oo) — 0. This equation is inte- 
grable and leads to an exact (energy-like) expression for 
the mid-point height 

^^^W{l^)^W{h), (6) 

where W{lTr) is the disjoining pressure for a wetting film. 
Making use of the Gibbs-Duhem equation dW {1.,^) / dh = 
Ztt in © , and reinserting the appropriate geometrical fac- 
tors for arbitrary tilt angles, we obtain 

, , S (seca-1) ^ , , 

fu,(/i,a)« \- sec a Dy, lT,(h) + , (7) 




FIG. 1: Schematic diagram of the wedge and cone geometries 
showing diverging lengthscales. 

where the ellipses denote vanishing contributions in the 
/i — > limit, and the critical amplitude takes the 
non-trivial value 

Dy, = ^- . (8) 

Thus, the next-to-leading-order correction is directly 
proportional to the planar wetting film thickness and 
becomes macroscopic in the same limit. The same ex- 
pansion Q) is also obtained in two dimensions from a 
transfer-matrix analysis and is valid in both the mean- 
field and the weak fluctuation regime. 

The divergent terms in the interfacial height {T)) con- 
tribute directly to the divergence of the excess adsorption 
in the wedge F. Thus, for vdW forces, F contains leading 
and next-to-leading-order contributions proportional to 
and h~^/^ respectively. Both these contributions to- 
gether with their corresponding critical amplitudes have 
been measured with remarkable accuracy in experimen- 
tal studies of Ar and Kr on silicon wedges by Mistura 
and co-workers 6j. 

It is also instructive to consider the scaling proper- 
ties of the equilibrium profile. Again, for convenience, 
we only consider the shallow wedge case, although the 
results can be readily generalised to more acute sys- 
tems. From the Euler-Lagrange equation, it follows that 
ri{x) — At„(x/^a;, x/^ii) where = ^a/h is roughly the 
point of contact of the meniscus with the wall. From this, 
one sees that the absolute interfacial height l{x) divides 
into two regions: a filled region for x < S^x corresponding 
to the central meniscus, and a tail for x > S^x where the 
relative interfacial height decays towards its planar value 
Ztt (see Fig. 2). The crossover between the regimes occurs 
over a distance comparable with the complete wetting 
transverse correlation length ^|| . Within the filled region, 
the meniscus is, to high accuracy, parabolic with 

Kx) ~ + ^ (9) 

whilst for X > S^x the profile decays exponentially quickly 
to the planar value 

l{x) « a\x\ + Z^(l-HCe~(^-«-)/«ii), (10) 
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with C a non-universal constant of order unity. These 
quantitative features wiU be of relevance below. 

We are now in position to study complete wetting in a 
cone geometry and establish our main result, the covari- 
ance with wedge filling for the mid-point heights. The 
substrate height is described by the function il}(x,y) = 
tana r, where r = ^ -V is the radial coordinate 
(see Fig. 1). As with the wedge, the cone strongly en- 
hances the adsorption compared to the planar wall and 
the divergence of as /i — > is accompanied by the 
divergence of the radial correlation length ~ and 
mid-point roughness ~ \J\vlIc. Thus, the transition 
is not fluctuation dominated and can be understood at 
mean-field level. In the cone, the macroscopic meniscus 
follows the surface of an inscribed sphere of radius 2S//i, 
since there are two equal radii of curvature contributing 
to the Laplace pressure. Therefore, the dominant contri- 
bution to the mid-point height is Ic ~ 2E(secQ; — 
twice the value of the corresponding wedge result. 

A more microscopic description is afforded by the in- 
terfacial Hamiltonian Q and its drum- head model gener- 
alisation. Minimisation of the Hamiltonian leads to the 
Euler-Lagrange equation for the relative height ?7(r) = 
l(r) — ar, 



(11) 



which is solved subject to the same boundary conditions 
?7'(0) = —a and r]'{oo) = 0. Unlike the much simpler 
wedge geometry, this equation is not integrable and the 
analogue of the energy-like equation ® reads 



dr 



(12) 



which is not an explicit expression for Ic- Despite the 
absence of integrability, one may still obtain the asymp- 
totic divergence of Ic analytically from consideration of 
the scaling properties of the interfacial profile. We are 
not aware that analytical results for this system have 
been presented before. From the Euler-Lagrange equa- 
tion, we have r;(r) = /7rAc(r/^r, ''/'?||), where = 2Y,a/h 
roughly separates the meniscus region for r < from 
the tail of the profile for r > ^r- Now, within the small 
angle approximation, the meniscus is very well described 
by 



l{r) 



Ic + 



4E 



(13) 



similar to the wedge result 0. The behaviour of the 
tails, however, is completely different for the cone. For 
distances r > , the profile decays algebraically slowly 

l{r) ^ ar+ (l - |) . (14) 

compared with the exponential decay for the wedge (|10|l . 
It follows that there exists no rescaling of the bulk or- 
dering field h which maps the full cone profile onto the 
wedge profile. 



Incidentally, the presence of the algebraic tails con- 
tribute significantly to the adsorption F in a cone, some- 
what analogous to critical adsorption phenomena. More 
specifically, in a cone of finite radius R, F contains a con- 
tribution linear in R arising from the tails in addition to 
the usual projected area term, that scales as R^. 

Now focus on the determination of the mid-point conic 
interfacial height Ic- The energy-like equation H12|) for the 
cone is different to that for the wedge © due to the pres- 
ence of the non-integrable final term, implying that both 
leading and next-to-leading order contributions to the di- 
vergence of lc{h) which are different to those in l^, eq. iQ. 
The leading order divergence differs due to a contribution 
— Ea^/2 in (|12|l from the integral over the meniscus re- 
gion r < ^r- Thus, the leading order divergence of Ic is 
twice that of Iw, consistent with the pure macroscopic 
considerations above. Similarly, the algebraic tails give 
rise to a contribution of order W{lTr) which alters the 
next-to-leading order singularity. Carefully making al- 
lowance for the full scaling properties of the profile and 
the specific form of the algebraic tail, it follows, after 
some algebra, that Ic diverges as 



lc{h,a) 



2E(seca-l) 



sec a Dc Kih) + 



(15) 



where again the ellipses denote vanishing contributions 
in the h —-^ limit. The next-to- leading-order correction 
term is once more proportional to the planar wetting film 
thickness but has a critical amplitude Dc given by 



Dc 



2^r 

1-/3,- 



(16) 



distinct from the wedge. It is only at this point that the 
covariance between the interfacial heights for the wedge 
and cone is manifest. Recalling that the planar wetting 
film thickness grows as h~^' , it follows that both diver- 
gent terms in Ic mimic precisely the wedge expression Q 
but at an effective bulk ordering field h/2. This leads 
directly to our central result eq. Numerical results 
for the mid-point heights obtained from minimisation of 
Hamiltonian @ in the wedge and cone geometries are 
shown in Fig. 2. It can be shown that covariance also 
applies to the full probability distribution functions 



Pc{l; h,a) =Pyj{l; - ,«) 



(17) 



where, in an obvious notation, Pc{l]h,a) denotes the 
PDF for the mid-point interfacial height in the cone, etc. 
By virtue of eq. and (|13|l , the covariance also extends 
to the shape of the interfacial profiles in the correspond- 
ing meniscus regions. All the above statements are valid 
for small values of h when the mid-point interfacial height 
is much larger than the planar film thickness. 

At this point, we make two remarks: 

• Note that the nature of the next-to-leading order 
correction term in the divergence of lw{h) and lc{h) re- 
sembles Deryagin's correction to the Kelvin equation for 
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h 



FIG. 2: Graphs of lw{h) (wedge), lc{h) (cone) and lc{2h) 
(shifted cone) obtained from numerical minimisation of the 
interfacial model for vdW forces. Dimensionless units, where 
A/X'' = Ea^ and h = h\/T,a^ are used. 

capillary condensation in a slit of width L the 
presence of thick complete wetting layers, condensation 
is shifted with respect to the bulk and occurs when the 
width satisfies 

L « I + 2DUh), (18) 
h 

where Deryagin's amplitude D takes the same value as 
Only for systems with short-ranged forces, for which 
(3^° = 0, do the right-hand sides of Q and ijTSJ) have 
a direct geometrical interpretation. The similarity is not 
coincidental and reflects disjoining pressure contributions 
to the capillary slit, wedge and cone free energies, respec- 
tively. 

• In addition to experimental studies, the present ex- 
ample of geometrical covariance can also be investigated 
in Ising model simulations. The reason for this is that, 
for systems with short-ranged forces, the expression H15|l 
for the divergence of the mid-point height also applies 
to pyramidal-like geometries where the height-function 
is 4^{x,y) = tana max(|a:|, \y\). Within the simple-cubic 
Ising model, both wedge [l^l and pyramidal-like geome- 



tries [13, with a = 7r/4, have been used previously to 
study the critical filling transition occuring at zero bulk 
field as 9 tt/A. Using the same geometries, but with 
stronger complete wetting surface fields (so that 6 — Q) 
and off bulk coexistence, one could readily test the ef- 
fective Hamiltonian prediction of covariance by directly 
comparing the simulation results for the two geometries. 

Finally, as stressed earlier, the covariance between the 
local interfacial heights at wedge and conic filling does 
not contradict the central theorem of density functional 
theory: the confining potential is uniquely determined by 
the density profile. What covariance implies, however, is 
that different confining potentials can lead to identical lo- 
cal interfacial properties albeit for rescaled fields. More 
generally, this implies that quite different examples of in- 
terfacial phase transitions share common superuniversal 
properties. This is perhaps more visible in the previously 
discussed covariance between 2D critical wedge filling and 
2D critical planar wetting. On approaching the critical 
filling phase boundary, — > a, at bulk coexistence, the 
mid-point wedge PDF satisfies Pw{l] 0, a) = Pti{1\ 6 — a), 
where PTr{l;0) is the planar critical wetting PDF writ- 
ten as a function of the contact angle. In these, and 
other examples of covariance, the relationships between 
the corresponding interfacial heights in the different ge- 
ometries only emerges after specific calculations involving 
both interfacial and more microscopic models. However, 
the fundamental physical origin of such relations in a 
more general context remains obscure. 

We believe that these different covariances point to- 
wards a more comprehensive theorem, perhaps within 
the framework of density functional theory, that would 
enable us to further understand the relation between ge- 
ometry and fluid interfaces. 
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